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ABSTRACT 

We present the results of the analysis of the internal velocity dispersions, a v , for the sample of 16 
distant galaxy clusters (0.17^3 z^, 0.55) provided by the Canadian Network for Observational Cosmology 
(CNOC). Different a v estimates are provided, all based on an interlopers removal algorithm, which is 
different from that originally applied by Carlberg et al. (1996). We find that all such methods provide 
a v estimates which are consistent within < 10% among themselves and with the original estimates 
provided by the CNOC collaboration. This result points in favor of a substantial robustness of currently 
applied methods for optical studies of the internal cluster dynamics. The resulting distribution of velocity 
dispersions is used to trace the redshift evolution of the cluster abundance with the aim of constraining 
the matter density parameter, f2,„. We find that constrains on fi TO are very sensitive to the adopted value 
of <7 8 = er 8 Sl^ (a — 0.4-0.5), as constrained by the local cluster abundance. We find that, as og varies 
from 0.5 to 0.6, the best fitting density parameter varies in the range 0.3^ f2 m S; 1.0. A further source 
of uncertainty in constraining O m is due to uncertainties in the correction for the a^-incompleteness of 
the CNOC sample. This calls for the need of better understanding the constraints from the local cluster 
abundance and increasing the statistics of distant clusters in order to suppress the systematics related 
to the sample completeness criteria. 

Subject headings: Cosmology: galaxies - clusters - theory - large-scale structure of the universe. 



1 INTRODUCTION 

Galaxy clusters represent the virialization stage of ex- 
ceptionally high peaks of initial density perturbations on 
comoving scales of about 10 /i _1 Mpc 7 . Therefore, their 
abundance is sensitive to the amplitude of the fluctuation 
power spectrum on such scales. Standard analytical meth- 
ods based on the approach originally devised by Press & 
Schechter (1974, PS hereafter) show that the number den- 
sity of clusters of a given mass provides a robust constraint 
on erg = CTgf^, where erg is the r.m.s. fluctuation ampli- 
tude within a sphere of 8/i -1 Mpc radius, Q m is the mat- 
ter density parameter and a ~ 0.4-0.5, weakly dependent 
on f2 m , on the presence of a cosmological constant term 
and on the shape of the power-spectrum (e.g., White, Ef- 
istathiou & Frenk 1993). 

However, while theoretical predictions provide the num- 
ber density of clusters as a function of their mass, obser- 
vations give the cluster abundance as a function of some 
observable quantity, like the A-ray luminosity, the A-ray 
temperature or the velocity dispersion of member galaxies, 
which are a posteriori connected to mass at different de- 
grees of reliability. Despite the variety of methods and of 
data sets employed to trace the local cluster abundance, 
all the analyses consistently indicate that <7g ~ 0.5-0.6 
(Eke et al. 1996; Kitayama & Suto 1997; Oukbir, Bartlett 



& Blanchard 1997; Girardi et al. 1998a; Pen 1998; Marke- 
vitch 1998; Borgani et al. 1999, B99 hereafter; Viana & 
Liddle 1999). Once a model is tuned so as to reproduce 
the abundance of local clusters, its evolution mainly de- 
pends on f2 m (e.g., Oukbir & Blanchard 1992). Therefore, 
having a statistical sample of high-redshift clusters with 
reliable mass determinations would allow in principle to 
break the degeneracy between Q m and erg. 

The growing availability of high-redshift cluster sam- 
ples selected in the A-ray band led in the last few years 
to a flurry of activity along this line. Using the redshift 
distribution of clusters from the Einstein Medium Sensi- 
tivity Survey (EMSS, Gioia et al. 1990), Sadat, Blanchard 
& Oukbir (1998) found consistency with a critical-density 
(f2 TO — 1) Universe (cf. also Reichart et al. 1999). B99 an- 
alyzed the ROSAT Deep Cluster Survey (RDCS, Rosati et 
al. 1995, 1998) and pointed out that current uncertainties 
in the evolution of the mass-luminosity relation prevent 
one from obtaining strong conclusions on f2 m even from a 
flux-limited sample as deep as RDCS (z^ 1). 

In this respect, the possibility of measuring A-ray clus- 
ter temperatures would circumvent this problem to a fair 
degree, since it is more directly connected to the mass than 
the A-ray luminosity. Indeed, Eke et al. (1998) analyzed 
the Henry (1997) sample of A-ray cluster temperatures, 
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extending out to z = 0.33. They concluded that the X-ray 
temperature function (XTF) favor a low-density Universe, 
with fl m ~ 0.4 ± 0.2. Quite remarkably, Viana & Liddlc 
(1999) analyzed the same data set and claimed that, af- 
ter accounting for systematics in both the data set and in 
the theoretical PS framework, a critical density Universe is 
still viable as far as the evolution of the XTF is concerned. 

It is clear that such ambiguities can in principle be elimi- 
nated by resorting to cluster data at a substantially higher 
rcdshift, where differences among different Fricdmann ge- 
ometries rapidly increase. However, the price to be paid 
in this case is that a much smaller number of clusters with 
highly accurate data is presently available (e.g., Luppino 
& Gioia 1995; Donahue et al. 1998; Rosati et al. 1999). 

An alternative way to trace the cluster mass function 
back in rcdshift is offered by measurements of the veloc- 
ity dispersions, a Vl of member galaxies. This represents a 
well established technique, which has been already exten- 
sively applied to samples of local clusters (e.g., Zabludoff, 
Huchra, & Geller 1990; Girardi et al. 1993; den Hartog 
& Katgcrt 1996; Fadda et al. 1996; Mazure et al. 1996; 
Borgani et al. 1997, and references therein). 

Systematic estimates of the velocity dispersions for a 
statistical sample of distant clusters has been performed 
for the first time by Carlberg et al. (1996, C96 hereafter). 
They analyzed the Canadian Network for Observational 
Cosmology (CNOC) sample of clusters (Yee, Ellingson & 
Carlberg 1996), which includes 15 X-ray selected clusters 
from EMSS and Abell 2390. Carlberg et al. (1997b, C97 
hereafter) used the CNOC sample to trace the evolution of 
the cluster abundance and concluded that fi m values be- 
low unity are in general to be preferred (see also Bahcall, 
Fan & Cen 1997). However, this result depends sensitively 
on (a) the z ~ normalization from the number density of 
local clusters with velocity dispersion above a given value, 
and (6) on the procedure to convert the distribution of 
cluster velocity dispersions into the distribution of clus- 
ter masses. In particular, one of the questions raised by 
C97 concerned the comparison between the different al- 
gorithms applied to estimate a v for local clusters and for 
the distant CNOC clusters. If the estimator used for local 
clusters provides a v which are biased upwards, this would 
give an overestimate of a$ and, therefore, of Q m . 

A reliable estimate of a v faces several problems, such 
as the presence of foreground and background interlop- 
ers, velocity anisotropy in galaxy orbits, the presence of 
substructures and the limited amount of data available. 
Girardi et al. (1993) showed that different methods pre- 
sented in the literature to estimate a v (Yahil & Vidal 1977; 
Zabludoff, Huchra & Geller 1990; Beers, Flynn & Geb- 
hardt 1990) give similar results on well sampled clusters, 
while only the robust estimator by Beers et al. (1990) 
seems to be efficient when only ~ 10 galaxy redshifts per 
cluster are available. Besides the estimate of u v from red- 
shifts of member galaxies, a further critical issue concerns 
the identification of genuine cluster members by removing 
foreground and background interlopers. 

For instance, Girardi et al. (1993) rejected interlopers 



by following a guess for the outer limits of the redshift 
range encompassed by each cluster. Two different, more 
refined, methods for interloper removal have been subse- 
quently developed by Fadda et al. (1996, F96 hereafter; see 
also Girardi et al. 1998b, G98 hereafter) and by den Har- 
tog & Katgert (1996; see also Adami et al. 1996; Mazure 
et al. 1996), both based on combining information on pro- 
jected position and velocity of each galaxy, although by 
using rather different procedures. Quite remarkably, such 
procedures of cluster member selection lead to velocity dis- 
persions which are in good agreement (cf. F96 and Adami 
et al. 1998). Significant differences were found only for 
a small number of clusters (8 out of 74) which were rec- 
ognized by the procedure of F96 as multiple structures in 
rcdshift space. 

Despite their differences, a common feature of all such 
methods is that they are based on the identification of 
individual interlopers to be removed. Quite differently, 
the method followed by C96 to estimate a v for CNOC 
clusters statistically subtract the mean density of field 
galaxies from the redshift space of the cluster. The result- 
ing velocity dispersions, which arc about 13% lower than 
without the background subtraction, have been shown by 
Mushotzky & Scharf (1997) to be in good agreement with 
those inferred from the X-ray temperatures provided. 
Furthermore, Lewis et al. (1999) have shown that virial 
masses for CNOC clusters are quite consistent with those 
obtained from X-ray data analysis. 

In this paper, we perform a complete re-analysis of the 
cluster velocity dispersion of the CNOC sample by ap- 
plying different algorithms which are similar to that used 
by F96 for local clusters. Such estimates of the velocity 
dispersions will be used to estimate the evolution of the 
cluster abundance through the redshift-dependence of the 
a v distribution. 

Therefore, the final goal of this paper will be to answer 
to the following two questions, (a) Does the explicit back- 
ground subtraction method, applied by C96 to estimate a v 
for CNOC clusters, provide consistent results to those from 
methods applied by F96 and G98 to local cluster samples? 
(b) Having a reliable determination of velocity dispersions 
for a sample with the statistics and the rcdshift extension 
of CNOC, are we able to provide robust constraints on f2 m 
from the evolution of the cluster abundance? 

The structure of the paper is as follows. In Section 2 we 
briefly describe the CNOC sample. In Section 3 we intro- 
duce the different methods to estimate a v . After showing 
the results of the application to the CNOC clusters, we 
compare the results to those provided by C96 (cf. also 
Carlberg et al. 1997a). In Section 4 we use such results on 
<j v to place constraints on er 8 and fl m . We draw our main 
conclusions in Section 5. 

2 THE CNOC SAMPLE 

A comprehensive description of the CNOC cluster sam- 
ple is provided by Yee, Ellingson & Carlberg (1996). In the 
following we only describe the main features which are rele- 
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vant to our analysis, while we refer to that paper for a more 
complete presentation. The CNOC cluster redshift survey 
includes all the 15 EMSS clusters (Gioia et al. 1990; Henry 
et al. 1992) which have L x > 4 x 10 44 erg s" 1 (h = 0.5 and 
q = 0.5) and / > fu m = 4x 10~ 13 ergs _1 cm -2 for X-ray 
luminosities and fluxes in the [0.3,3.5] keV energy band, 
redshift in the range [0.17,0.55] and -15° < 5 < +55° for 
the declination. The solid angle covered by EMSS in this 
declination range is of 587 sq. deg. Also included in the 
original CNOC sample is Abell 2390, which however does 
not belong to EMSS. Although it is included in our a v 
analysis, we will not use it for the purpose of constraining 
cosmological models. 

The CNOC surveys provided a fairly large number (50 to 
200) of galaxy redshifts per cluster over a region covering 
up to 6/i _1 Mpc projected on the sky. Galaxy magnitude 
selection criteria have been chosen so as to minimize con- 
tamination due to foreground/background objects. This 
makes the CNOC clusters very well suited to follow the 
redshift evolution of the cluster internal dynamics. 
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Fig. 1. — The velocity-space galaxy density, as provided by the 
adaptive-kernel reconstruction method. Units on the x-axis 
are velocities in 10 3 kms _1 , while the y axis is in arbitrary 
units. 



3 Analysis of the CNOC sample 

In order to select member galaxies and compute veloc- 
ity dispersions, a v , we apply the same procedure adopted 
by F96 and G98 for nearby clusters. We remind that the 
other recent a v analysis for an extended sample of nearby 
clusters (Mazure et al. 1996) leads to velocity dispersions 



which are in good agreement, although it resorted to a dif- 
ferent procedure for the selection of cluster members (cf. 
F96 and Adami et al. 1998). 

3.1 Cluster Member Selection 

The member selection for the CNOC clusters is per- 
formed by considering the whole sampled region around 
each cluster, with the exception of MS1512+36, which we 
analyze only within lft, _1 Mpc. This cluster appears as 
very elongated and shows strong complexity in the two 
dimensional galaxy distribution. Therefore, we prefer to 
avoid distant cluster regions where this complexity makes 
difficult a dynamical analysis (see below). The initial list 
of candidate cluster members do not include galaxies which 
are very far from the cluster redshift range and, therefore, 
are clearly foreground/background objects. 

The identification of cluster members proceeds in two 
steps. 

Firstly, we perform the cluster membership selection in 
velocity space by using only redshift information. We use 
the adaptive kernel method (Pisani 1993) to find the sig- 
nificant (> 99% c.l.) peaks in the the velocity distribu- 
tion. Only galaxies belonging to such peaks are consid- 
ered as candidate cluster members. In Figure 1 we plot the 
velocity-space galaxy density for all the CNOC clusters, as 
provided by the adaptive-kernel reconstruction method. 

From the above member selection analysis only the 
MS0906+11 cluster, shows two significant peaks, which 
are separated by > 3000 km s -1 in velocity space (here- 
after MS0906+lla and MS0906+llb). According to our 
procedure, these two peaks show a significant overlapping 
(> 40%) and they are not perfectly separable, thus lead- 
ing to a cluster with uncertain dynamics. The amount of 
evidence of this binary nature (e.g. C96; Carlberg et al. 
1997a), as well as the strong disagreement between a v of 
the whole system (1723 km s -1 from the present analysis 
and 1893kms _1 from C96) and the A-ray temperature 
reported by Lewis et al. (1999; 8 keV, corresponding to 
~ 1150 km s -1 under the hypothesis of energy equiparti- 
tion between gas and galaxies), indicates that this cluster 
should be better considered as composed of two separated 
peaks. 

Also the velocity distribution of MS1512+36 is very 
complex. The cluster peak is surrounded by secondary 
peaks, which have a smaller density and are separated from 
the main peak. This complexity further motivates our 
choice to estimate a v for this cluster only within 1 /i _1 Mpc. 

All those galaxies which pass through the first velocity- 
space selection are analyzed in the second step, which uses 
the combination of position and velocity information. This 
procedure, which is based on the "shifting gapper" algo- 
rithm, proceeds in two iterative steps. Firstly, galaxies are 
assigned to radial bins of 0.4 fe _1 Mpc width or larger, in or- 
der to contain at least 15 objects. Secondly, those galaxies 
that, within each bin, are separated in cluster rest-frame 
velocity by > 1000 km s -1 from the main body of the ve- 
locity distribution are identified as interlopers. These two 
steps are iterated until the number of cluster members 
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does not change anymore. Figure 2 shows the plots of rest- 
frame velocity versus projected cluster-centric distance for 
those clusters where the "shifting gapper" found interlop- 
ers (indicated with the open circles). We note that the 
MS1358+62 cluster contains a close system, which corre- 
sponds to a southern group identified in velocity space and 
already noted by C96 (the vertical line mark the separa- 
tion between these two systems). In order to exclude this 
secondary structure, hereafter we analyze the MS1358+62 
cluster only within 1.2 /i _1 Mpc. 
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Fig. 2. — Rest-frame velocity versus projected cluster-centric 
distance for each cluster where the "shifting gapper" succeeded 
in rejecting at least one galaxy. Distances on the x-axis are in 
units of /i _1 Mpc and the rest-frame velocities on the j/-axis are 
in units of 10 3 kms -1 . Interlopers are indicated with the open 
circles. As for the MS1358+62 cluster, the thin vertical line 
marks the separation between the main cluster and a nearby 
subgroup (see text). 

Finally, we estimate the corrections for velocity gradi- 
ents, which are however always quite small (e.g. Girardi 
et al. 1996). We verify that only A2390 has significant 
(at > 99% c.l.) velocity gradients. For this cluster we ap- 
ply a correction by subtracting the velocity gradient from 
each galaxy velocity and renormalizing the velocities so as 
to have their mean velocity unchanged. This correction 
results in a decrease of a v (see below) by 58kms _1 . 

3.2 The velocity dispersion estimation 

We compute robust estimates of (line-of-sight) veloc- 
ity dispersions, cr ro 6, by using the ROSTAT routines by 
Beers, Flynn, & Gcbhardt (1990), after applying the rel- 
ativistic correction and the usual correction for velocity 



errors (Danese, De Zotti, & di Tullio 1980). We use the 
bi-weight scale estimator when more than fifteen mem- 
ber galaxies are available, while the gapper estimator is 
applied otherwise (see also Girardi et al. 1993), 

Several authors (F96; den Hartog & Katgert 1996; Gi- 
rardi et al. 1996) have shown that velocity dispersion pro- 
files for individual clusters are often a strong, either in- 
creasing or decreasing, function of the radius in the central 
cluster regions. This behavior may be due both to velocity 
anisotropics and to the dark matter distribution, although 
it is not easy to disentangle between these two effects (e.g., 
Merritt 1987). However, most integrated profiles, a v (< R) 
(i.e., a v evaluated by using the all the galaxies within the 
projected radius R), become flat in the external cluster 
regions (F96; Girardi et al. 1996). This suggests that pos- 
sible velocity anisotropies does not significantly affect the 
value of o~ v when it is estimated over a wide cluster region, 
so as to be representative of the total kinetic energy. 

Figure 3 shows the integrated velocity dispersion pro- 
files, <r v (< R), for all the CNOC clusters. The clusters 
which are sampled up to external cluster regions really 
shows a flattening of the profiles (e.g., A2390, MS0451+02, 
MS1621+26), while the behavior of profiles for some others 
is still uncertain (e.g., MS1006+12). 
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Fig. 3. — Integrated line-of-sight velocity dispersion profiles, 
where the dispersion at a given radius is a rob estimated by 
considering all galaxies within that radius. The bootstrap er- 
ror bands at the 68% c.l. are shown. Distances on the a;-axis 
are in units of fe _1 Mpc, and the velocity dispersions on the y- 
axis are in units of 10 3 kms -1 . As for MS0906+11 the solid and 
dotted lines correspond to peak a and b, respectively (errorbars 
are omitted here for reasons of clarity). 
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In order to test the robustness of the a v computations, 
we also use three other estimators: (a) a weighted ver- 
sion of the robust estimate <r w -rob, (b) the usual standard 
deviation er s , and (c) its weighted version a w - s - As for 
the weighted estimators, to each galaxy is associated a 
magnitude-dependent geometric weight, which has been 
introduced in order to account for the geometry of the 
field surveyed around each cluster and the magnitude com- 
pleteness of the sample (cf. Yee, Ellingson, & Carlberg 
1996 for further details). The result of the a v analysis 
are reported in Table 1 where we give for each cluster: 
the total number of galaxy redshifts, N to t, available for 
each cluster (Column 2); the number of member galax- 
ies found in peaks, N m , as recognized by our interloper- 
removal algorithm (Column 3); the bi-weight mean cluster 
velocity (Column 4); the values of a ro b with the relative 
bootstrap error (at 68% c.L; Column 5); and the other 
three estimates of the velocity dispersion that we consid- 
ered (Columns 6-8). 




Fig. 4. — The comparison between line-of-sight velocity dis- 
persions (in units of kms" 1 ) as estimated by Carlberg et al. 
(1996), <tc96, and the four types of estimates we obtain in this 
work (see text). The la bootstrap errorbars are shown only 
in the <JC96-a T ob comparison. The two very discrepant points 
refer to the two sub-clusters in which MS0906+11 has been 
divided. 

In Figure 4 we compare our estimates of a v to those, 
CTC96, provided by C96 and based on the explicit back- 
ground subtraction algorithm. Apart from the exception 
of MS0906+11, which we divide into two subclumps, the 
present a v estimates show an overall agreement with ac96 ■ 
We find that the median values of the ratio a v /o'C96 for the 



different estimators lie in the range 0.96 — 1.04. In partic- 
ular, it is cr ro b/(7c96 = 1-04I_q;q7, where errors are at 90% 
c.l. It is remarkable that different methods give a v which, 
on average, agree within ^ 10%, thus confirming the re- 
sult by Girardi et al. (1993) that different estimators of a v 
give statistically similar results on cluster samples. This 
supports the reliability of such cluster dynamical studies 
from optical observations. 

Unless otherwise specified, we adopt in the following the 
a ro b estimator, which is the most directly comparable to 
that used by F96 and G98. 

We also show in Figure 5 the correlation between a v and 
bolomctric X-ray luminosity for CNOC clusters (filled cir- 
cles). This correlation is also compared to that obtained 
for local clusters by combining the XBACs sample (Ebel- 
ing et al. 1996) with the G98 sample of velocity dispersions 
(open circles). 




Fig. 5. — The Lboi-& v relation for local (open circles) and dis- 
tant CNOC clusters (filled circles). For the local data, we show 
results for those clusters of the Girardi et al. (1998) sample, 
which are also included in XBACs (Ebeling et al. 1996). We re- 
stricted to those clusters whose a v is estimated at least with 30 
galaxy redshifts. For CNOC clusters, X-ray luminosities have 
been taken from the EMSS (Henry et al. 1992). Short-dashed, 
long-dashed and continuous lines are the direct, inverse and 
bisector linear regression for the local data points. The light 
horizontal and vertical dashed lines indicate the Lx,boi and a v 
limits above which we select CNOC clusters for deriving cos- 
mological constraints. 

Bolomctric luminosities for CNOC (XBACs) clusters, 
Lx,boi, have been estimated by multiplying luminosities 
in the [0.3-3.5] keV ([0.1-2.4] keV) energy band by a 
temperature-dependent bolometric correction factor. This 
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TABLE 1 
Cluster Velocity Dispersions 



Name 


Ntot 


N m 


V 


°Vo6 


"w — rob 




&w — s 








kms -1 


kms -1 


kms -1 


kms 


kms" 1 


A2390 


252 


191 


68608 


12621 H 


1213 


1306 


1024 


MS0016+16 


65 


41 


164729 


H27±lf| 


1042 


1029 


1014 


MS0302+16 


58 


27 


127333 


not 1 * 


688 


678 


630 


MS0440+02 


78 


47 


58793 


715^ 


815 


709 


542 
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factor has been computed under the assumptions of pure 
bremsstrahlung ICM emission and power-law approxima- 
tion for the Gaunt factor, g(E,kT) oc (E/kT)-°- 3 , which 
is quite accurate for kT > 2 keV (cf. B99). Tempera- 
tures are available for all but two CNOC clusters, namely 
MS1231+15 and MS1621+26 (Lewis etal. 1999). For such 
two clusters they have been estimated from the velocity 
dispersion according to kT — (cr„/350kms _1 ) keV. As for 
the double cluster MS0906+11, we assigned the same over- 
all temperature, kT = 8 keV, to both subclusters. The re- 
sulting bolometric correction factors range from 1.6 to 3.7. 
According to Fig. 5, it turns out that there is no appre- 
ciable evolution for the Lx.boi~o~v relation, thus consistent 
with the lack of significant evolution in the Lx,boi~T rela- 
tion at z^ 0.4 (e.g., Mushotzky & Scharf 1997) and also in 
agreement with the recent analysis by . The only outlier is 
MS1512+36, which is one of the two highly substructured 
objects that we discussed above. The solid line represents 
the bisector between the direct and the inverse log-log lin- 
ear regression (short- and long-dashed lines, respectively) 
for local clusters. The resulting Lx~<J v relation is given 
by 

(1) 



logLx = 5.1t^log<7 w -14.2 



H3.0 



with a scatter A^ x — AlogLx/logLv — 0.36 around 
the best-fitting relation and in agreement with the recent 
analysis by Wu, Xue & Fang (1999). The upper and lower 
errors represent here the difference with respect to the di- 
rect and the inverse linear regression, respectively. Eq.(|l|) 
has been obtained by considering the 53 clusters of the 



G98 sample, which have a v > 500 km s -1 and at least 
30 galaxy members (cf. G98). If Tx oc a 2 , as expected 
under the assumption of isothermal gas and hydrostatic 
equilibrium, then eq.(Q) would imply Lx.toi oc T£ 2 5 for 
the luminosity-temperature relation, a results which is not 
far from recent calibrations of the Lx,boi~Tx relation for 
rich clusters (e.g., Allen & Fabian 1998; Arnaud & Evrard 
1999, and references therein). 

4 Constraining Vt m 
4-1 Modeling the <r v distribution 

The starting point of this analysis is represented by the 
Press-Schechter (1974, PS hereafter) expression for the 
comoving number density, at redshift z, of virialized halos 
with mass in the range [M, M + dM] : 



N(M, z) dM 



P S c {z) 

IT M 2 



exp 



0~M 

Sc(z) 2 



d log a M 



d log M 
dM . 



(2) 



Here p is the present-day average matter density and 
5 c {z) is the linear-theory overdensity extrapolated at the 
present time for a uniform spherical fluctuation collaps- 
ing at redshift z. It is convenient to express it as S c (z) = 
8 C fi(z)[D(0)/D(z)], where D(z) is the linear fluctuation 
growth factor (see, e.g., Peebles 1993). For a critical- 
density Universe, <5 c .o = 1.686 with a weak dependence 
on fi m . We take for Q m < 1 the expression provided by 
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Lacey & Cole (1993) and Kitayama & Suto (1996). The 
quantity om is the present day linear r.m.s. fluctuation 
within a sphere of radius R = (SM/Anp) 1 / 3 and is speci- 
fied by the power-spectrum, P(k), of density fluctuations. 
We assume for P(k) the parametric CDM-like expression 
provided by Bardeen et al. (1986). Accordingly, its profile 
and amplitude are determined by the shape-parameter F 
(~ 0.2 from the the observed galaxy distribution; e.g., Efs- 
tathiou, Bond & White 1992) and by erg, the r.m.s. fluctu- 
ation amplitude within a sphere of 8/i _1 Mpc radius. The 
reliability of eq.(||) for predicting the mass distribution of 
virialized halos has been tested and debated at length in 
the literature (e.g., Lacey & Cole 1994; Bryan & Norman 
1998; Gross et al. 1998; B99; Governato et al. 1999) and 
we refer to such papers for further details. 

In order to convert the distribution of cluster masses 
into that of one-dimensional velocity dispersions, a v , we 
use the relation 



M(a v ) = 



1129 km s- 1 



A c (z) 
178 



-1/2 



© > 



(3) 



where E(z) = [(l + z) 3 n m + {l + z) 2 (l-Q m -Q A ) + n A ]^ 2 
and A c is the ratio between the average density within 
a virialized halo and the critical background density at 
z (e.g., Bryan & Norman 1998). It is A c = 178 for 
Q m = 1 while the expression for £l m < 1 has been taken 
from Kitayama & Suto (1996). Eq.(||) corresponds to 
the case of an isothermal halo density profile for f a = 1 
(Binney & Tremaine 1987), with slightly smaller values, 
0.9^ f a ^S 1, for steeper profiles at the virial radius (e.g., 
Navarro, Frenk & White 1996). Since different halos have 
different dynamical histories, we expect that an intrinsic 
scatter, Am = AM/M, should be present around eq.(||). 
Once the M-a v relation and its scatter are calibrated, the 
distribution of velocity dispersions, N(a v ), is obtainable 
by convolving eq.(^) with a Gaussian distribution having 
dispersion A M . 

In order to estimate f ffv and Am, we resort to TV-body 
simulations, based on the Adaptive P3M code by Couch- 
man (1991). We analyze two sets of simulations, which 
have been run for an Einstein-de-Sitter (EdS) model and 
for a flat low-density model with fL„ — 0.4 (AO. 4). Both 
models have been chosen to be consistent with the abun- 
dance of local clusters. Each simulation box is 250 /i _1 Mpc 
aside and five realizations have been run for each model. 
Such simulations have already been used by B99 to test 
the PS mass function and we refer to that paper for fur- 
ther details. Once clusters are identified, f a and Am are 
estimated from their M-a v virial relation. The results 
of this analysis are reported in Figure 6, where we plot, 
at three different redshifts, f„ and Am as a function of 
the cluster number density, lower n c ; corresponding to the 
population of richer clusters. We find that f a and Am 
are essentially independent of both the simulated models 
and the evolutionary stage. Only the EdS result on f a 
at z — 0.2 for the smallest n c i (i.e., for the most mas- 



sive cluster population) shows a marginal departure from 
the general behavior. We checked that this is generated by 
one single cluster, appearing in one of the five realizations, 
which has an anomalous high <j v . Based on this result, we 
adopt in the following the best fitting values, f a = 0.93 
and Am = 0.15. 
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Fig. 6. — The normalization / CT (upper panels) and the scatter 
Am (lower panels) for the M-a v relation [cf. eq.^] from N- 
body simulations are plotted against the cluster number den- 
sity. Left and right panels are for the EdS and the AO. 4 model, 
respectively (see text). In each panel, circles, squares and tri- 
angles refer to the outputs at z = 0, 0.2, 0.6, respectively. Er- 
rorbars, which are reported only for the z = 0.2 output, are la 
Poisson scatter. 



4-2 Maximum-likelihood analysis 

In order to constrain cosmological parameters we com- 
pare the cluster distribution on the (a v ,z) plane with 
model predictions for samples having the same complete- 
ness and selection criteria of the CNOC sample. A simi- 
lar approach has been recently used by Eke et al. (1998) 
to follow the evolution of the cluster X-ray temperature 
function for the Henry (1997) sample, and by Borgani et 
al. (1998) to follow the evolution of the X-ray luminosity 
function in the ROSAT Deep Cluster Survey (Rosati et al. 
1998). 

This method is based on partitioning the (a v ,z) plane 
into elements of size da v dz and on computing the model 
probability \(a V: z)da v dz of observing a CNOC cluster 
with velocity dispersion a v at redshift z, given the com- 
pleteness and selection criteria of the sample. If the bin 
width is small enough that such probabilities arc always 
much smaller than unity, then the likelihood function C 
of the observed cluster velocity dispersions and redshift is 
defined as the product of the probabilities of observing ex- 
actly one cluster in dz da v at each of the (cr V) i, Zi) positions 
occupied by the CNOC clusters, and of the probabilities 
of observing zero clusters in all the other differential ele- 
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ments of the (a Vl z) plane. Assuming Poisson statistics for 
such probabilities, we obtain 



£ = Y[ [H a v,i, z%, )dzda v ( 

i 

-X(a V; j ,Zj)dzda v 



— A(<7 V f i ,Zi)dz da v 



x 1> 



(4) 



where the indices i and j run over the occupied and empty 
elements of the (a v ,z) plane, respectively. Defining, as 
usual, S = — 21n£, and dropping all the terms which do 
not depend on the parameters of the cosmological model, 
we find 

S = -21n£ = -2j2HH<rv,i,Zi)} 

i 

2 / dz / da v X(cr v , z) , (5) 



where the integral is performed over the redshift and ve- 
locity dispersion intervals spanned by the cluster sample 
(see also Marshall et al. 1983). 
The function X(a v , z) is given by 



X(a v , z) = N(a v , z) fsky , 



(6) 



where N(a v , z) is the comoving number density of clusters 
with velocity dispersion a v at redshift z, dV(z) is the co- 
moving volume element in the redshift interval [z, z + dz] 
and f s ky is the sky-coverage, i.e. the effective fraction of 
the sky surveyed by EMSS. In order to account for boot- 
strap errors in the a v estimates we compute N(a v ,z) by 
convolving the Press-Schechter-based predictions with a 
10% Gaussian-distributed relative scatter, which is repre- 
sentative of the statistical uncertainties in a v (cf. Table 1). 
Since f s k y for EMSS clusters depends on the detect-cell 
flux, wc compute it for given a v and z by integrating over 
all the possible fluxes, thereby accounting for the scatter 
in the Lb i-a v correlation of eq.(Q). In order to account 
for the fact that detect-cell fluxes are smaller than total 
cluster fluxes, we divide the latter by a factor 2.1, which 
has been shown by Henry et al. (1992) to be appropriate 
at least for clusters in the redshift range [0.22-0.33] (cf. 
their Figure 1). 

A further selection effect that we have to account for in 
the estimate of the likelihood function is connected to the 
fact that, at a given a v , the X-ray luminosity- and flux- 
limit criteria cause faint clusters to be missed in the CNOC 
sample. In order to correct for this a\,-incompleteness, we 
follow a similar procedure to that adopted by C97. We se- 
lect those CNOC clusters which have a v > 800 km s -1 and 
Lx,boi > 10 45 ergs _1 . Then we resort to a complete local 
cluster sample to estimate the fraction J-goo of all clusters 
with a v > 800 km s -1 , which also satisfy the above lumi- 
nosity selection. Finally, we correct for the a v incomplete- 
ness by multiplying by J^^qq the sum over the occupied 
cells in the r.h.s. of eq.(||). 

To estimate ^-"sooj we cross-correlate the ENACS cluster 
sample (Katgert et al. 1998), which at z^S 0.1 is represen- 
tative of the whole cluster population for a v > 800 km s -1 



(Mazure et al. 1996), with the XBACs (Ebeling et al. 
1996), which includes all the Abell/ACO clusters identi- 
fied in the ROSAT all-sky survey and is 80% complete for 
fluxes /p. 1-2.4] > 5 x 10~ 12 ergs _1 cm -2 . Since at z = 0.1 
this flux corresponds to a luminosity L x , [0.1-2.4] — 2 x 
10 44 ergs -1 , we are reasonably guaranteed that, for typical 
bolometric corrections, XBACs contains all the Abell-type 
clusters with L x ,bol > 10 45 ergs _1 in the redshift range 
covered by ENACS. Therefore, Tsoo is given by the frac- 
tion of ENACS clusters with a v > 800 km s -1 and belong- 
ing to the XBACs with Lx,boi > 10 45 ergs _1 . It turns out 
that ENACS contains 26 clusters with a v > 800 km s -1 , 
out of which 6 have Lx,boi > 10 45 ergs _1 . Therefore, we 
obtain = 4.3 ± 2.0, where the uncertainty is the la 
Poissonian error. We note that using local data to estimate 
J-goo is consistent with the lack of significant evolution in 
the Lx.boi-&v relation shown in Fig. 5. 

In their analysis, C97 used the data on Lx and a v 
by Edge & Stewart (1991) and found = 2.0 ± 1.0. 

The consequence of a smaller J-%J is that a smaller frac- 
tion of clusters is expected to be lost below the luminos- 
ity limit. Therefore a smaller number of clusters with 
a v > 800 km s -1 is expected at high redshift and a larger 
value of fi m would be implied. This will be discussed on 
a more quantitative ground in §4.3 below. 

Best estimates of the model parameters are obtained 
by minimizing S and confidence regions are estimated by 
allowing for standard increments AS. Unless otherwise 
specified, in the following we will restrict our analysis only 
to the case of vanishing cosmological constant. 

4-- 3 Results 

Figure 7 shows constraints on the a^-^m plane taking 
r = 0.2 for the power-spectrum shape (we checked that 
results are weakly dependent on T within the range al- 
lowed by galaxy clustering data; e.g., Peacock & Dodds 
1994). Such results are based on adopting the unweighted 
robust estimator for a v (cf. Section ^ and taking only 
clusters with a v > 800 kms -1 and Lx,boi > 10 45 ergs _1 . 
The plotted iso-likclihood contours corresponds to la, 2a 
and 3cr confidence levels. The dashed curves are the a§- 
£l m relation coming from the local cluster abundance and 
correspond to a$ = 0.50, 0.55 and 0.60, while the shape 
of the (78-f2 m relation is that provided by Girardi et al. 
(1998a). 

Both the small number of clusters on which the anal- 
ysis is based and the limited leverage in redshift allow a 
rather high level of degeneracy in the a$-Q m plane. As 
expected, more stringent constraints on Q m come by com- 
bining the results from CNOC clusters with those from 
the local cluster abundance. Taking cfg = 0.55, wc find 
fl m = 0.65^q Jg (at the 2a c.l. for one significant param- 
eter), thus ruling out a critical density model at a high 
confidence level (if fl\ = 1 — fl m , then this result changes 
into Q m = O.53IQ12). However, the variation of the best- 
fitting f2 OT as ag is increased from 0.5 to 0.6 is much larger 



8 



than internal statistical uncertainties; indeed, we find that 
n m = 0.35±g;o? and 1.05±g;?2 for a s = 0.5 and 0.6, respec- 
tively. Therefore, although uncertainties from the local 
cluster abundance are believed to be rather small, they 
propagate into large uncertainties in the determination of 
tt m from distant cluster data (cf. also Colafrancesco, Maz- 
zotta & Vittorio 1997). This is due to the fact that CNOC 
clusters are rather rich systems and, therefore, probe the 
high-mass tail of the mass function, which is very sen- 
sitive to the power spectrum normalization. On the one 
hand, this confirms the good news that results from a small 
number of clusters with precise determinations of a v can 
be used to constrain Q m , once they are combined with lo- 
cal data. On the other hand, the bad news is that current 
uncertainties in the local cluster abundance are still large 
enough to prevent placing strong constraints on the den- 
sity parameter. We regard this as a general problem one 
has to face with, even assuming that the analysis of high z 
sample is perfectly under control. 



0.8 



0.6 - 



0.4 



0.2 




7. — Iso-likelihood contours on the (ag,il m ) plane for 
0.2. Results refer to the reference analysis method, 
Iso-likelihood contours are for 



Fig. 

r = 

as described in Section 
AS = 2.30,6.17 and 11.8, which corresponds to la, 2a and 
3a c.l. for two significant fitting parameters. The three dashed 
curves are the a$-£l m relation from the local cluster abundance, 
for as = 0.50, 0.55 and 0.60, from lower to upper curves, while 
the scaling is taken from Girardi et al. (1998b). 

A further source of uncertainty in constraining Q m arises 
from the correction for the a v -incompleteness. We find 
that, as jFgoo varies, within its la range, from 2.3 to 6.3 
(cf. §4.2), the best-fitting Q, m decreases from ~ 0.9 to 
~ 0.5. Therefore, the effect of this uncertainty is also non- 



negligible and calls for the need of a substantially larger 
number of clusters with both Lx and a v determinations 
to suppress the Poissonian uncertainty in the estimate of 

Furthermore, we note from Fig. 5 that there are two 
clusters, MS1224+20 and MS0906+lla, whose values of 
L>x,boi are only slightly below the adopted luminosity limit. 
Owing to residual uncertainties in the bolometric correc- 
tion, one may wander by how much results would change if 
such clusters were included in the analysis. In this case, we 
find that, for <7 8 = 0.55, the best-fitting f2 m only decreases 
by 0.06. 

In order to clarify the crucial role of the low-z nor- 
malization for tracing the redshift evolution of the cluster 
abundance, we compare in Figure 8 the comoving cluster 
number density, with a v > 800kms _1 , to model predic- 
tions for a critical and an open low-density model, assum- 
ing different normalizations from the local cluster abun- 
dance. The estimate of the CNOC cluster abundance is 
obtained by dividing the nine selected clusters into two 
redshift intervals, 0.15 < z < 0.30 and 0.30 < z < 0.55, 
which contain five and four clusters, respectively. Num- 
ber densities are then estimated by applying the 1/V max 
method (Avni & Bahcall 1978) and following the same 
procedure as C97. Apart from the details of the model 
constraints, which are in general agreement with the re- 
sults from Fig. 7, this plot gives a visual impression of 
why a change in the normalization from the local cluster 
abundance turns into a change of the preferred £l m value. 
For instance, n(> 800, z = 0) ~ 3 x 10~ 6 ( /i^Mpc)" 3 , 
as coming from the virial analysis of local clusters (e.g., 
Mazure et al. 1996; Fadda et al. 1996), would be consis- 
tent with a critical density Universe, while n(> 800, z = 
0) ~ 10 _6 ( /i -1 Mpc)~ 3 , as implied by the estimates of the 
local A -ray temperature function by Eke et al. (1996) 
and Markevitch (1998), would instead favor a low-density 
Universe. 

Since this is not the first analysis aimed at constraining 
Q m and as from the CNOC cluster abundance evolution 
(C97; Bahcall et al. 1997), it is worth comparing our re- 
sults to those from such previous analyses. C97 provided 
the number density of clusters with mass within a physi- 
cal radius of 1.5/i Mpc above a given limit by dividing 
the CNOC sample into two redshift bins. It is interest- 
ing to note that, despite the value of J-goo used in that 
analysis is about twice as large than ours, their results 
prefer lower 17 m (cf. Bahcall et al. 1998, for a stronger 
claim for a low density parameter, Q m ~ 0.3, from CNOC 
clusters). The reason for this difference is mainly due to 
the fact that C97 use a fixed normalization for the lo- 
cal cluster abundance. In particular, they noticed that, 
if ENACS velocity dispersions are thought to be overes- 
timated by 13%, then the implied cluster abundance is 
compatible with that found by Eke et al. (1996) from the 
A -ray temperature function. In fact, the number den- 
sity of local clusters with a v > 800 km s -1 , as given by 
C97, would be compatible with ag = 0.5. In view of the 
good agreement we found here between velocity disper- 
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Fig. 8. — The redshift evolution of the comoving number density of CNOC clusters with a v > 800kms _1 , compared with model 
predictions. Left and right panels refer to a critical-density and to an open model with fi m = 0.4, respectively. In each panel, the 
three different curves refer to three different normalizations from the local cluster abundance. We assume F = 0.2 for the shape 
parameter of the power spectrum. The filled dots are the CNOC cluster number density in two redshift bins (see text). Vertical 
errorbars are la Poissonian uncertainties in n(> a v ), while horizontal errorbars indicate the extent of the redshift bins. 



sions estimated with the explicit background subtraction 
and with the method applied to ENACS clusters, we argue 
that there is no evidence that ENACS velocity dispersions 
should be overestimated. 

5 Conclusions 

We analyzed the internal velocity dispersions, a v , of 
CNOC clusters (e.g., Yee, Ellingson & Carlberg 1996), by 
applying the algorithm originally developed by Fadda et 
al. (1996, F96; cf. also Girardi et al. 1998b, G98) for 
the analysis of local clusters. After removing interlopers, 
we applied four different a v estimators. By using the ro- 
bust estimator, we found that for a v > 800 km s -1 and 
Lx,boi > 10 45 ergs -1 the CNOC sample contains nine clus- 
ters, on which we base our analysis. In order to account for 
the clusters that we miss with respect to an ideal sample 
complete for a v > 800 km s -1 , we resort to local data in 
order to reliably estimate the fraction J-soo of clusters that 
in such a sample have Lx,boi > 10 45 ergs -1 . The resulting 
a v are used to trace the evolution of the cluster abun- 
dance in the 0.17^S z^S 0.55 redshift range probed by the 
CNOC sample, with the purpose of placing constraints on 
the a%-Q, m plane. Such constraints can then be combined 
with results from the local cluster abundance to constrain 
the cosmological density parameter. 

The main results of our analysis can be summarized as 
follows. 

(a) The explicit background subtraction method applied 
by Carlberg et al. (1996, C96) on CNOC clusters 



provides a v estimates which are fully consistent with 
those provided by the method applied in this paper. 
For instance, we find that the median of the ratio 
between the robust estimator and the C96 results is 

&rob/PC96 = I-U4_ 07 . 

(b) Current uncertainties in the a^-Slm relation from the 
local cluster abundance are still large enough not to 
allow CNOC data to distinguish to a high confidence 
level among a low-density Universe with n m ~ 0.3 
and a critical density scenario. 

Far from meaning that the evolution of the cluster abun- 
dance can hardly be used to place significant constraints on 
cosmological parameters, the results of our analysis any- 
way indicate that the situation is more complex than some- 
times suggested. Somewhat surprisingly, the availabil- 
ity of precise cluster mass measurements at high redshift 
through both a v and X-ray temperature measurements 
may be only of partial help in increasing the strength of 
the constraints unless (a) systematic uncertainties in the 
local mass-function are reduced at a 10% level, and (b) 
the completeness criteria of high-z cluster samples are well 
under control. 

We acknowledge stimulating discussions with Neta Bah- 
call, Andrea Biviano, Alain Blanchard, Isabella Gioia, 
Michael Gross, Simon Morris and Joel Primack. We thank 
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body code. SB and MG acknowledge SISSA in Trieste for 
its hospitality during the preparation of this work. 
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